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ABSTRACT 

PSR J2021+3651 is a 17 kyr old rotation powered pulsar detected in the radio, 
X-rays, and 7 -rays. R powers a torus-like pulsar wind nebula with jets, dubbed the 
Dragonfly, which is very similar to that of the Vela pulsar. The Dragonfly is likely 
associated with the extended TeV source VER J2019-I-368 and extended radio emission. 

We conducted first deep optical observations with the GTC in the Sloan r' band to search 
for optical counterparts of the pulsar and its nebula. No counterparts were detected 
down to r' > 27.2 and > 24.8 for the point-like pulsar and the compact X-ray nebula, 
respectively. We also reanalyzed Chandra archival X-ray data taking into account an 
interstellar extinction - distance relation, constructed by us for the Dragonfly line of 
sight using the red-clump stars as standard candles. This allowed us to constrain the 
distance to the pulsar, D = 1 . 8 l }4 kpc at 90% confidence. R is much smaller than the 
dispersion measure distance of ~12 kpc but compatible with a 7 -ray “pseudo-distance” 
of 1 kpc. Based on that and the optical upper limits, we conclude that PSR J2021-I-3651, 
similar to the Vela pulsar, is a very inefficient nonthermal emitter in the optical and 
X-rays, while its 7 -ray efficiency is consistent with an average efficiency for 7 -pulsars 
of similar age. Our optical flux upper limit for the pulsar is consistent with the long- 
wavelength extrapolation of its X-ray spectrum while the nebula flux upper limit does 
not constrain the respective extrapolation. 
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Subject headings: stars: neutron - pulsars: general - pulsars: individual: PSR 

J2021+3651 


1. Introduction 


The 104 ms pulsar J2021+3651 was discovered in the radio with the Arecibo telescope in a 
deep search for radio pu lsations towards unid entified ASCA X-ray sources spatially coinciding with 
EGRET 7 -ray objects (jRoberts et al.l 120021 1. With the characteristic age Tc ~ 17 kyr and spin- 
down luminosity E ~ 3.4 x 10^® erg s“^ this pulsar is among the youngest and most energetic 
rotation-powered pulsars known. The canonical dipole magnetic field estimated from the pulsar 
period and its derivative is R ~ 3.2 x 10^^ G. A tent ative detection of 7 -ra y puls ations with 


the pulsar period in the EGRET data was reported bv iMcLaughlin fc CordesI (j2003). Later 7 - 


ray o bservations with AGILE satellite (jHalpern et al.ll2008l l and Eermi observatory (jAbdo et al 


2009|) hrmly established a double-peaked pulse p rofile and a power -law spectrum with a photon 


index T ~ 1.7 and cutoff energy of ~ 2.9 GeV (jAbdo et al 


20131 1. PSR J2021-k3651 was also 


identified in X-rays with Gha ndra, and weak pulsations, also with t he double-peaked profile, were 
detected at a 4a signihcance (jHessels et al.l 120041 : lAbdo et al.l [2009|) . The pulsar X-ray spectrum 
contains thermal and n qn-thermal components from the surface and magnetosphere of a neutron 
star (NS), respectively (jVan Etten et al.ll2008l ). Ghandra also revealed an extended pulsar wind 
nebula (PWN) G75.2-1-0.1 whose brightest internal part, within ~ 30^' of the pulsar, has a torus¬ 
like mor phology with axial jets . By its specihc spatial shape this PWN was dubbed the Dragonfly 
Nebula (jVan Etten et al.ll2008l l. A fainter diffuse emission is extended up to several arcminutes. 


The PSR J2021-I-3651 position is projected on the Gygnus-X region, one of the richest known 
regions of star formation in the Galaxy. A bright extended TeV source MGRO J2019-I-37 was 
identified wit h the Milaqro sky survey in this region with a 20 TeV flux of 80% of those of the 
Crab Nebula (jAbdo et al.l 120071 ). The source was suggested to be associated with the Dragonfly, 
which was recently confirmed by observations with the VERITAS observatory. VERITAS resolved 
the source into two objects (jAliu et al.ll2014lb The brightest one, VER J2019-I-368, has a hard 
spectrum resembling the spectrum of Vela X - a TeV PWN system powered by the Vela pulsar. 
VER J2019-I-368 coincides also with an extended region of non-thermal radio emission. 

The most controversial parameter of P SR J2021-h3651 is the distance. The NE2001 model for 
the Galactic distribution of free electrons (jCordes fc Laziol 120021 ) for the pulsar line of sight (/ = 
75?21. h = 0?11) and the dispersion measure DM ~ 370 pc cm“^ yield a distance D ~ 12 kpc (e.g., 
Roberts et al.l 120021 ). Comparing a hydrogen absorbing column density obtained fro m first X-ray 


observations and the total Galactic HI column density along the pulsar line of sight, iHessels et al 


^Based on observations made with the Gran Telescopio Canarias (GTC), instated in the Spanish Observatorio del 
Roque de los Muchachos of the Institute de Astrofsica de Ganarias, in the island of La Palma, programme GTC3-11B. 
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(|2004l) suggested D Ri 10 kpc. I Van Etten et alJ (j2008|) performed similar analysis of subsequent 
deeper X-ray observations and found the distance of 3- 4 kpc. The pulsar polarization rotation 
measure implies a minimal ~ 5 kpc (|Abdo et al.ll2009l b Adopting the latter value a nd assuming 
that PSR J2021-I-3651 was born near the center of VER J2019-I-368, Aliu et al. ( 2014 1 estimated a 
possible transverse velocity of th e pulsar to be ^ 8 40 km s“^, which is about 3 times higher than 
the average for pulsar velocities (jHobbs et al.ll200f)l ). Finally, the distance can be as low as 1.5 kpc 
if the pulsar is located within the Cyg nus-X region. This locatio n is consistent with the empirical 
7 -ray “pseudo-distance” relation le.g.. lSaz Parkinson et al.ll2010l l suggesting Z1 ~ 1 kpc. 


By many multiwavelength properties PSR J2021-I-3651 and its PWN are similar to the Vela 
pulsar plus PWN system, but, in contrast to the Vela, this pulsar has never been studied in the 
optical. We report first deep optical observations of the PSR J2021-I-3651 field performed with the 
10.4-m Gran Telescopio Canarias (GTC). We also address the issue of the distance discrepancies, 
using the Chandra archival X-ray data and red-clump stars as standard candles, and compare the 
optical results with the X-ray ones. The details of observations and data reduction are described 
in Sect. [21 our results are presented in Sect. Eland are discussed in Sect. IH 


2. GTC data 

2.1. Observations and data redaction 

The pulsar field was imaged with the Optical System for Imaging and low-intermediate Res¬ 
olution Integrated Spectroscopy (OSIRI^) in the Sloan r' band at the GTG on 2011 September 
28. Sixteen dithered 158-second exposures were obtained using the OSIRIS standard image scale 
of Cf.'254/pixel with the field of view of 7f8 x 7f8. The field was exposed on a mosaic of two GGDs, 
with the target source placed on CCD2. The observations were carried out during dark time, and 
the conditions were photometric, with seeing varying from Cf.'8 to l''.l. 

Standard data reduction and analysis, including bias subtraction and flat-fielding, was per¬ 
formed with IRAF tools. To eliminate shifts between individual exposures, we collected a set of 
unsaturated stars in the field and aligned the images to the one with the best seeing using IRAF 
routines. The alignment uncertainty was <0.1 pixel. All exposures were then combined and yielded 
a final image with a mean seeing of Cf.'9, airmass of 1.8, and total integration time of ~ 2.5 ks. 


^For instrument details see http://www.gtc.iac.es/instruments/osiris/ 
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Fig. 1.— Top: ~ 60" X 60" GTC/OSIRIS 
optical Sloan r' (top-left) and Chandra/ACIS- 
S X-ray (top-right) image fragments of the 
PSR J2021-I-3651 field. The pulsar is marked 
by “-b” in the X-ray image and contours indicat¬ 
ing the X-ray PWN boundary and pulsar loca¬ 
tion region are shown in both images. Bottom: 
13" X 13" fragment of the Sloan r'-band de¬ 
picting the pulsar vicinity. The circle shows 3a 
X-ray pulsar position uncertainty. Two sources 
nearest to the pulsar are labeled. 


2.2. Astrometric referencing and photometric calibration 

For astrometric referencing, positions of 10 suitable astrometric standards from the USNO- 
Bl.O catalogu^ were used. Their pixel coordinates were measured on the combined image with 
the IRAF task imcenter. The IRAF ccmap routine was applied to the astrometric transformation 
of the image. Formal rms uncertainties of the astrometric fit for the combined image were ARA < 
[f!l7 and ADec < C/.G3. Accounting for the nominal catalog uncertainty of ~ ff!2, this results in 
conservative estimates of la referencing uncertainties of < C/.'26 for RA and < (y.'24 for Dec. 

The photometric calibration was carried out with G158-100 Sloan standard observed the same 


®http://www. nofs.navy.mil/data/fchpix/ 
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night as our target. The atmospheric extinction of 0.10 ±0.01 mag airmass ^ for the Sloan r'-band 
was taken from the OSIRIS user manual. The resulting magnitude zero-point for our r' image is 


29±13±0±02. 


3. Results 


3.1. Searching for the Dragonfly and pulsar optical counterparts 


The 60" X 60" pulsar field fragment of the GTC/OSIRIS r'-band image, which contains the 
brightest part of the Dragonfly Nebula, is shown in the top-left panel of Figure [H The image 
is smoothed with a Gaussian kernel with width a = 2 pix. It is compared with the respective 
Chandra/KCISS 0.5-8 keV X-ray image (top-right panel of Figured]), obtained by merging all 
the available archival datc0 (112 ks effective exposure in total). The data were reprocessed with 
the CIAO v.4.6 chandra-repro tool with GALDB v.4.5.9. The X-ray image is binned by two ACIS 
pixels, smoothed with one-pixel Gaussian kernel, and shown in log-intensity scale. The X-ray PWN 
is comprised of a SW jet, NE counter-jet, and two arcs, which are oriented perpendicular to the 
jets. The arcs are bel i eved to be associated with the PWN equatorial torus seen almost edge 
on (jHessels et ahl 120041 : 1 Van Etten et ahl l2008l j . Contours in the X-ray image indicate the outer 
boundary of the torus-like PWN, where it blends with the background, and the region around the 
pulsar (marked by the cross). In the top left panel of Figure [1] the X-ray contours are overlaid 
on the optical image. The vertical bold strip crossing the left side of the optical image is a bleed 
line from a bright over-saturated background star located outside the fragment. Two wave-shaped 
horizontal curves near the top side of the image are detector artifacts. Comparing the optical and 
X-ray images, we do not find any significant extended optical feature correlated with the X-ray 
morphology of the compact Dragonfly PWN. However, at an arcminute scale comparable to the 
faint diffuse X-ray emission extent, there are some background variations containing bright and dark 
regions seen in all individual exposures as well . Examination of Ha images of the field from the 
INT Photometric Survey of the Northern Sky (iBarentsen et ahl 120141 ) shows that these variations 
correlate with the Ha emission variations. 

The immediate pulsar vicinity is enlarged in the left-bottom panel of Figure [TJ where the 
r ^2^4 circle is cent ered at the pulsar X-ray position with RA = 20:21:05.46 and Dec = ±36:51:04.8 
(jHessels et ahl 120041 1. It corresponds to the 3 cj pulsar position uncertainty which accounts for the 
optical astrometric referencing and pulsar X-ray position uncertainties. No significant point-like 
objects are detected within the pulsar error circle. The closest reliably detected point-like source 
“a” with r' = 24.40 ± 0.04 is located at about 4"8 or at ~ 6cr from the pulsar X-ray position. More 
distant object “b” has a lower brightness of 25P01 ± 0P05 and is located at ~ l/.M or about 12cr 


^PI Roberts, C/iarirfra/ACIS-S, OBsID 3901, date of obs. 2003-02-12, 19 ks exposure; PI Romani, Chandra/ACIS- 
S, OBsIDs 7603 and 8502, date of obs. 2006-12-29 59-1-34 ks exposure 
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from the pulsar. Because of their large offsets, both sources are unrelated to the pulsar. 

Using our optical image we, therefore, can set only upper limits on the pulsar and the Dragonfly 
Nebula flux densities in the Sloan r' band. For the pulsar, we used a mean background deviation 
within a circular aperture with a 4 pixel (~ 1^') radius centered at the pulsar position. We accounted 
for an aperture correction of OUl derived using bright background stars. The resulting 3a upper 
limit on the pulsar flux density is ^0.04 /rJy (r' >27.2). For the nebula, we used an elliptical 
aperture with semi-axes of 6"2 and IC/.T and a position angle of 137° centered at the pulsar, which 
encapsulates most of the X-ray PWN equatorial torus emission. The 3a upper limit on the spatially 
integrated flux density of the PWN is <0.36 /iJy (r' > 24.8). 


3.2. Distance and interstellar extinction 


It is possible to construct an extinction-distance relation for the direction towards the pul- 
sar utilizing red-clump stars as st andard candles, follow i ng a method described, for instance, in 
Lopez-Corredoira et all (120021 ) and ICabrera-Lavers et al.l (j2005l ) . The method was used previously 
to constr ain distances and e xtinctions for several sources. Some examples are the X-ray binary 4U 
1608—52 (iGiiver et al.ll2010ll. six anomalou s X-ray pulsars (IDurant &: van Kerkwiikll2006l l. and two 
7 -ray pulsars (jPanilenko et al.ll2012l . l2013l l. 


In the top panel of Figure O we show K vs J — K band color-magnitude diagram for stars 
from the 2MASS All-Sky Point Source Catalogu^ located within 0?3 of the pulsar position. The 
red-clump (RC) branch, as well as main-sequence (MS) and asymptotic-giant (AGB) branches, are 
indicated. We divided the diagram into several magnitude bins, and in each bin we fitted J—K color 
distribution with a mixture of two Gaussians corresponding to the MS and RG branches. The AGB 
stars were eliminated by omitting all points located right of a boundary starting at J — K = 1.7 for 
small magnitudes and ending at 2U5 for large ones. The derived J — K colors of RG stars with their 
unce rtainties were then transf ormed into distances and interstellar extinctions Ay, using relations 
from iRieke &: Lebofskvl (119851 1 , assuming the absolute magnitude and the intrinsic color of the red- 
clump stars to be Mk = —1.62 ± 0.03, and {J—K)q = 0.68 ± 0.07, respectively (see the above-cited 
papers for details). Extinctions were transformed to hydrogen absorbin g column densities Ny usi ng 


a standard empirical relation Xh = Ay x (1.79 ± 0.03) x 10^^ cm ^ (jPredehl X: Schmittlll99,5l ) 


The resulting IVh- distance dependence is shown in the bottom panel of Figure [2j Nn increases 
with the distance reaching a limit of (15±4) x 10^^ cm“^ at distances ^ 5 kpc. Within uncertainties 



®see http://irsa.ipac.caltech.edu/applications/DataTag/, DataTag = ADS/IRSA.Gator7/2014/0814/100517_7624 
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J-K 



Fig. 2.— Top: K vs J — K diagram for the stars from the 2MASS All-Sky Point Source Catalogue 
located within 0?3 of the pulsar position {I = 75?22, b = 0?11). The main-sequence (MS), red- 
clump (RC), and asymptotic-giant (AGB) branches are indicated. Solid line shows a smoothed 
spline approximation to the RC stars mean colors; dashed segments are its extrapolations to higher 
and lower K magnitudes. Light-shaded region bounded by dot-dashed lines contains 95% (2cj) of 
RC stars. Bottom: Empirical A^H“distance relation for the PSR J2021-I-3651 direction derived using 
the RC stars colors from the diagram at the top panel. Bars are la uncertainties. The solid line and 
gray filled region are smoothing spline approximations to the data points and their uncertainties, 
respectively. They are linearly extrapolated to small and large distances (dashed lines). 
















is considered as not reliable at b < 5°. 


3.3. X-ray spectral analysis 


We reanalyzed the archival Chandra data in light of the A^H-distance relation. We extracted 
the pulsar spectra from all three C/iandra/ACIS-S sets using an aperture with a radius of C/.'Td 
centered at the pulsar position applying the CIAD v.4.6 specextract tool. We also extracted the 
PWN spectrum from an elliptical region with semi-axes of 6"2 and lOi'G and a position angle of 
137°, that encloses most of the PWN equatorial torus emission. The circle aperture of a 2" radius 
around the pulsar was excluded from this region. Backgrounds were taken from regions free of any 
sources on the ACIS-S3 chip, where the Dragonfly was exposed in all three Chandra data sets with 
live times of 19, 59, and 43 ks. Total numbers of source counts are ~ 1270 for the pulsar and ~ 
5250 for the PWN. 


To evaluate likelihoods, we used the statistics. To model the pulsar spectrum, we applied 
an absorbed sum of the power-law (PL) and thermal components. Any single component did not 
describe the data. For the thermal co mponent, we tried blackbody (BB ) and magneti c neutron 
star hydrogen atmosphere models NSA ( Pavlov et al. 1995 1 and NSMAX ( Ho et al. 2008 1. For the 
interstellar abs orption, we used the XSPEC photoelectric absorption phabs mode l with default cross- 
sectio ns bcmc (|Balucinska-Church &: McCammonlll992l l and abundances angr (jAnders &: Grevesse 
1989l l. 


To model the contribution of the PWN to the spectrum extracted from the pulsar aperture, we 
added second PL to the pulsar spectral model and fitted the PSR and PWN spectra simultaneously 
in the 0.3-10 keV spectral range. The second PL component photon index was tied with the PWN 
photon index, and Nh was set as a global parameter. Doing this, we also took into account the ratio 
of the PWN flux within the pulsar aperture to the total PWN flux of ~ 0.05, as it was estimated by 
Van Etten et al.l (j2008l l and independently confirmed by us via modeling of Chandra/ACIS PSF. 


The Ah “distance relation and its uncertainty were approximated by smoothing splines, shown 
by the line and gray filled region in the bottom panel of Figure El This relation wa s then used as 


a Ba.v esian prior information for the subsequent spectral fitting procedure (see, e.g.. lGelman et al 


20031 . for details). Technically, we assumed that for each distance the Ah value follows a Gaussian 
distribution with the mean and a taken from the approximations. We then run Markov chain Monte 
Carlo (MCMG) using the Go odman-Weare algorithm implemented as a python package emcee by 
Foreman-Mackev et al.l (j2013l ). For each model we kept 1000 steps after initial burn-in, whi ch is 


large enough considering that typical autocorrelation time (see, e.g. iGoodman fc WeardEOlflU was 


of order of several tens (50-90) of iterations. As 100 MCMG walkers ([Goodman fc Weard l2010l ) 
was used, we obtained 100000 samples in total. 


Posterior median values of spectral parameters with 90% credible intervals for the BB-I-PL and 
NSMAX-|-PL models are presented in Table [TJ The goodness-of-fit test (y^ values are in Tabled]) 
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shows that both models are equally consistent with the data. We present fit results for only one 
of hydrogen atmo sphere models, NSMAX 1260, which corresponds to the surface magnetic field 
B = 4 X 10^^ G ( Ho et al.l 120081 ') . The thermal component can be equally well fitted by NSA or 
by any other model from the NSMAX set; the resulting parameters do not depend appreciably on 
the choice of specific atmosphere model. We prefer more modern NSMAX models because they 
account for the partial ionization in stellar atmospheres, especially important at low temperatures 
(as in our case), while the NSA models were constructed for fully ionized NS atmospheres. 


Our results are generally co nsistent with those prese nted bv lHessels et al.l (j2004l l and IVan Etten et al. 


(|2008l b However, in contrast to IVan Etten et al.l (|2008l b we did not fix Nh at a best-fit value ob¬ 


tained from a separate analysis of the PWN spectrum. Nevertheless, the resulting IVh is d efined 


mainly by the PWN spectrum, and consistent with the one obtained by I Van Etten et al.l (120081 ) 


within uncertainties. Thereby it weakly depends on a particular model used to describe the pulsar 
spectrum and on a particular N^^-D relation. The distance D = 1.8^|' 4 kp<i§ is now mainly deter¬ 
mined by Vh and the adopted N-n-D relation. In this approach, only two parameters are defined 
by the thermal component: the radius R and the effective temperature T of the emitting area. 
Importantly, we were able to infer the absolute value of the radius R, not only the R/D ratio as 
it would have been without accounting for the Vh — D relation. As it is typical for pulsars where 


Pavlov et al. 

2001; 

Kirichenko et al. 

2014) 


factor of 2.5 larger than those for the hydrogen atmosphere model. For the latter, R = 12^^g km 
implies that emission can come from the bulk of the NS surface with the effective su rface temper¬ 


ature , redshifted for distant observer, T = 63^8 eV, close to that of the Vela pulsar ( Pavlov et al 
20011 ). For the BB model, R = I.SIJ'q km, which is compa tible with a canonical pulsar hot polar 
cap radius of ~ 0.6 km for a 100-ms pulsar (jSturrocklll97ll ). 


'Here and below we discuss the largest D range from Table [T] 





























Table 1: Posterior medians for the pulsar and the PWN equatorial torus region spectra. The BB+PL and NSMAX+PL models 
are for the pulsar spectrum, while the PWN is described by the PL model f- 


Model 

Nn 

10^^ cm“^ 

r 

-L psr 

^psr ^ 

10“^ photons 
keV“^ cm“^ s“^ 

T 

eV 

R 

km 

D 

kpc 

r 

pwn 

^pwn 1 

10“® photons 
keV“^ cm“^ s“^ 

(dof) 

BB+PL 

cr q+0.5 
^•0-0.5 

1 q+0.6 

-*-•'“'-0.6 


-1 0+1.5 
-‘-•'J-1.0 


-*-■^-0.1 

q 7+1.0 
y-'-0.9 

459(469) 

NSMAX+PL 

0-U_o,5 

1 

U.D_o 4 


12 0+^^-^ 
-'-^•'“'-9.6 


1 4+0-1 
-1-+-0.1 

10 

l'“'-'“'-0.9 

472(469) 


is the absorbing column density. The temperatures T and emitting area radii R for the BB and NSMAX spectral 
components are given as measured by a distant observer. For the NSMAX component, the gravitation redshift 
I + z = {1 — 2.952M/_R)“° ®, where M and R are the NS mass and circumferential radius in the Solar mass and km 
units, respectively, is fixed at 1.21. This corresponds to a canonical NS with M — 1.4 and it = 13 km. amj 

Tpsr.pwn are PL normalizations and photon spectral indexes for the pulsar (psr) and PWN (pwn), respectively. All 
errors correspond to 90% credible intervals derived via MCMC. 

^Pulsar fluxes in 2-10 keV range are logP^®’^ = —IS.S^Q g, —IS.djiJ'j [ erg cm“^ s“^ ], for BB+PL and NSMAX+PL 
models, respectively. PWN flux in the same range is logP^"“ = —12.2lg'j [ erg cm“^ s“^ ] for both models. 
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3.4. Multiwavelength spectra of the pulsar and PWN 


The best-fit A^h values obtained from the X-ray spectral analysis suggest a total interstellar 
extinction towards the Dragonfly Ay ~ 3.3 in the V band, whi ch results in the exti nction A^' ~ 2.8 
in the r' band using a standard extinction law with Ry = 3.1 (jCardelli et al.lll989l b Based on this, 
upper limits on the dereddened flux densities for the pulsar and PWN in the r' band are about 0.57 
/iJy and 4.85 //Jy, respectively. In Figure [3] we compare these limits with unabsorbed X-ray spectra 
of the pulsar (top panel) and PWN (bottom panel), fitted by BB-I-PL and PL models, respectively. 
For the PWN, the optical and X-ray data were obtained from the same spatial region enclosing its 
X-ray equatorial torus emission (see Sect. 13.1113.31 and Figure [I]). The solid line in the top panel of 
Figure [3] shows the total best-fit model, including the contribution of PWN nonthermal photons to 
the spectrum extracted from the pulsar aperture. The dashed line shows solely the PL component 
of the pulsar. As seen, the PWN contribution is substantial only in the high-energy tail. 

As seen from Figure [Sj the pulsar optical flux upper limit does not exceed the extrapola¬ 
tion of the best-fit X-ray spectral model to the optical. This i s typical for rotation powered pul 


2010 


sars detected in the optica l and X-rays (jPanilenko et al.l 120111 : iDurant et al.l l201ll : iMignani et al 


Shibanov et al.ll2006l L For all of them, the nonthermal component dominates in the optical. 


However, this component usually shows a break between the optical and X-rays with a signifi¬ 
cant spectral flattening in the optical. Our data do not exclude the presence of such a break for 
PSR J2021-I-3651, although the extrapolation of the X-ray PL component is still rather uncertain 
and the optical limit is not deep enough. The NSMAX-I-PL model, which equally well fits the X-ray 
data, does not change these conclusions. 

All torus-like PWNe, which have been detected in b oth spectral domains , also show spectral 
flattening in the optical in comparison with X-rays (e.g., IZharikov et al.ll2013l i. However, for the 
Dragonfly the situation is currently even less certain than for its pulsar. The nebula optical flux up¬ 
per limit overshoots the low-energy extrapolation of its X-ray spectrum (bottom panel of Figure [3]), 
and the presence of the break in the spectrum of this PWN remains an open question. 


_ The dereddened upper limits were obtained using the Ay — Ah rela tion of iPredehl fc Schmitt 

( I 995 I ). There exist other empirical relations of that kind. For instance, Gfiver fc Ozell ( 2nf)dl gi^ 


Ah = Ay X (2.21 ± 0.09) x 10^^ cm“^, also consistent with the results of iGorensteinI ( 1975l L 
Using this relation instead, we get smaller Ay ~ 2.7 and hence dereddened upper limits a smaller 
by a factor of 1.6, which does not change general conclusions of this Section. 


4. Dicussion 

Our rather deep, down to r' ~ 27.2, GTC optical imaging of the Dragonfly Nebula field 
allowed us to set upper limits on the optical flux densities of PSR J2021-I-3651 and its PWN. The 
non-detection of this energetic system can be attributed to high interstellar extinction towards 
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Fig. 3.— Top: Unabsorbed spectrum of PSR J2021+3651. The solid line is the best-fit model 
for the Chandra X-ray spectrum which includes BB and PL pulsar spectral components and the 
PWN contribution to the spectral extraction aperture (see text for details). The best-fit model 
is extrapolated towards the optical. The GTC dereddened 3a flux upper limit in the r' band is 
shown by the bar with the arrow. Dash-dotted and dashed lines with light- and dark-gray regions 
are the BB and PL pulsar spectral components with their 90% uncertainties, respectively. The 
difference between the solid and dashed lines is clearly visible at the high-energy tail and reflects 
the PWN contribution. Bottom: The solid line and gray region are the best-fit PL model of the 
X-ray spectrum of the PWN equatorial torus region with its 90% uncertainties, respectively. The 
GTC dereddened optical 3a flux upper limit for this region in the r' band is also shown. The 
error-bar crosses in each panel are the unfolded Chandra data. 
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the object, which is roughly about of 3™ in the r' band. Considering the data for other pulsars 
detected in the optical and X-rays, we conclude that 1-2 magnitude deeper optical observations are 
necessary to detect this pulsar and to reveal the expected spectral break between the optical and 
X-rays. This is, in principle, feasible with 8-10 m ground-based telescopes, such as the GTC, using 
a few hour exposure in r' at good seeing conditions, although observations at longer wavelengths, 
less affected by the interstellar absorption, would be preferable. Such observations would also be 
useful to better constrain the optical-X-ray spectral properties of the PWN. 


Our Xn-distance relation, constructed using the red-clump star method and compiled with the 
X-ray spectral analysis, supports previous suggestions that the pulsar is likely to be substantially 
closer to us than it is inferred from DM and the NE2001 model of the Galactic distribution of 
free electrons. Our est imate D = 1.8 j ’ .] kpc is compatible, within uncertainties, with the 3-4 
kpc range suggested bv IVan Etten et al.l (120081 ). However, our allowed distance range is shifted to 
lower distances. It suggests the association of the pulsar with the Gygnus-X region, located within 
2 kpc from the Sun, and is consistent with the y-ray “pseudo-distance” of ~ 1 kpc provided by 
the Fermi data. The reduced distance we found makes feasible parallax and pulsar proper motion 
measurements with VLBI. A possible source of systematic errors in our distance determination 
method originates from ambiguity in Ay — Ah relations, as stat ed above. Reprocessing the analysis 
of Sections [32] and [331 using the relation o uGiiver fc Ozel[ (|2009li . we obtain an even smaller distance 
of 1.3lj f kpc. 

Gomparing the DM of 370 pc cm“^, or electron column density Ag ~ 1.14 x 10^^ cm~^, with 
Ah of 6x10^^ cm“^ leads to an average ionization ratio of 19% along the pulsar line of s ight which 
is not too much larger than the 10% ionization found on average (e.g., [He et al.[[2013[ . figure 1). 
On the other hand, the NE2001 electron density model in the pulsar direction gives much smaller 


Ne = 0.7li:2x 


1020 

cm ^ for distance range of D = 1.8l^'4 kpc. There exist several indications 


that the NE2001 model s trong l y underestimates A^, in the vicinity of the PSR J2021-I-3651 direc¬ 
tion ( Gamilo et al. 2009 . 


201 ^ 


Arumugasamv et al.l [20141 b For instance, there is another pulsar 
J2022-I-3842 at only 1?8 from PSR J2021-I-3651 with very high DM = 429 pc cm~^. for whic h 
NE2001 model gives obviously overestimated distance D > 50 kpc (j Arumugasamv et al.l l2014l ). 
This may imply that there ar e dense clouds in t he Gygnus-X region, which are not taken into 
account in the NE2001 model ( Roberts et al. 2002). 


As for the X-ray thermal emission component of PSR J2021-I-3651, we cannot state definitely 
whether the thermal emission comes from a hot polar cap or the bulk of the NS surface. This 
depends on which, BB or a hydrogen atmosphere, model is applied to describe the thermal emission. 
The phase-resolved spectroscopy would be useful to distinguish between the two possibilities. We 
note, however, that if the thermal emission originates from the entire surface of the star (the case 
of the atmospheric model), PSR J2021-I-3651 has a rather small surface temperature for its age. 
According to the NS cooling theories, such a small temperature can not be explained by a standard 
cooling scenario, but it can be reached if the effects of superfluidity in the stellar interiors are 
invoked. This is possible, for instance, if the powerful direct Urea process of neutrino emission 
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Fig. 4.— X-ray and optical luminosities, Lx and Lopt, and respective efficiencies, r/x and rjopt, 
for pulsars of different charac t eristic age r detected in both spectral domains. The data are 
adopted from iDanilenko et al.l ([201j). Different pulsars are marked by different symbols. The 


PSR J2021-I-3651 data, derived in this work, are included (marked by asterix). 


operates in the star but is suppressed by superfluidity (|Yakovlev &: PethickI l2004l l . In another 
interpretation the direct Urea process is not allowed, but the cooling is enhanced due to the specific 
process of the neutrino emission accompanying the Cooper pair formation i n the neutron triplet 
superfiuid. Thi s is th e essence of the so-called minimal cooling scenario (j Gusakov et al.l l2004l : 


Page et al.ll2004l . 1200911 . 


The pulsar’s 0.1-100 GeV 7 -ray luminosit y ^ 5.9x 10^^ erg s ^ and efficiency rj.y = L^/E ~ 
1.8, derived i n the 2 nd Fermi Pulsar Catalog (jAbdo et al.ll201,3l ) using the distance of 10 kpc from 


Hessels et al.l (j2004l ) appear to be unreasonably high, and place the pulsar at the highest end of L.y 
and distributions of 7 -ray pulsars. In contrast, for the distance D = 1 . 8 lJ '4 kpc inferred from 
our analysis, logL.y = 35 . 315*3 [ ] and log 7 / 3 , = — 1 . 2 l'j *'3 become consistent with the average 

values of the respec tive distributions for 7 -ray pulsars with similar E and/or characteristic age (cf.. 
Figures 9 and 10 of Abdo et al. 2013 b 


The pulsar’s unabsorbed nonthermal X-ray flux in 2-10 keV range, derived from the X-ray 
spectral fits, is logFx = —13.5lQg [ erg cm“^ s“^ ]. The respective X-ray luminosity and efficiency 
are log Lx = 31.ll^^ [ ] and logr^x = ~5.4llg, assuming the distance range derived in this 

work. 


For the 90% distance upper limit D = 3.5 kpc, upper limits on the optical luminosity and effi- 
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ciency in the V band, assuming a flat spectrum, are logLopt ^ 29.9 [ erg s“^ ] and logr/opt ^ —6.7, 
respectively. In Figured] we compare the obtained X-r ay and optical efficienc ies and luminosities 
with the data for other pulsars observed in both ranges (|Danilenko et al.ll2013l l. According to Fig¬ 
ured! we can conclude that PSR J2021-I-3651, like the Vela pulsar, is inefficient in these ranges as 
compared to other substantially younger and older pulsars. We thus obtain a new member of the 
small sample of Vela-like pulsars forming a puzzling minimum in the optic al and X-ray efficienc y 


dependences on age at a characteristic age of ~ 10 kyr noticed pr eviously bvIZharikov et al.l ([2006j). 


No such minimum is visible in the respective y-r ay dependence (lAbdo et al.ll2013l l. Together with 
strong glitches and high polarization in the radio ( Hessels et al. 2004h. a bright NS t hermal emission 
component in X-r ays, a double arc X-ray PWN with lets (IVan Etten et al.ll2008l l. 7 -ray activity 
(lAbdo et al.ll2009l l. and association with a TeV source (lAliu et al.ll2014l l. this makes the Dragonfly 
pulsar and PWN remarkably similar to the Vela pulsar and its PWN. 
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